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An approximate analytic solution is constructed for the 2-d space charge limited emission by a 
cathode surrounded by non emitting conducting ledges of width A. An essentially exact solution 
(via conformal mapping) of the electrostatic problem in vacuum is matched to the solution of a 
linearized problem in the space charge region whose boundaries are sharp due to the presence of a 
strong magnetic field. The current density growth in a narrow interval near the edges of the cathode 
depends strongly on A. We obtain an empirical formula for the total current as a function of A 
which extends to more general cathode geometries. 
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The study of space charge limited (SCL) current, ini- 
tiated in the beginning of the last century [1,2], contin- 
ues to be of great current interest [3-9]. These works 
are important for the design of high power diodes, tech- 
niques of charged particles beams, physics of non-neutral 
plasmas including plasma sheath, and other numerous 
applications. The modelling of SCL ionic flow in cylin- 
drical and spherical geometry [3] is also necessary for the 
incrtial-electrostatic confinement of fusion plasmas. Un- 
fortunately only the planar 1-d case permits an analytic 
solution [1,2] and as pointed out in a recent review [5] 
even "the seeming simple problem of 2-d planar emission 
remains unresolved" . This motivated the present work 
which provides a semi-analytical solution for a prototype 
2-d model similar to that studied in [6]. We obtain for 
the first time, we believe, a reasonable analytic approx- 
imation for the currents at the edge of the cathode - an 
important (though usually undesirable) feature of SCL 
diodes [6,7]. An extension of our methods should facili- 
tate dealing with this problem to any desirable accuracy 
thus providing an alternative to PIC simulations. 

Model. We consider the current between two conduct- 
ing flat electrodes where the anode, whose potential is V, 
is an inflnite plane separated by a distance D from the 
grounded cathode which is an infinitely long strip paral- 
lel to the anode. Our assumptions are: 1) The cathode 
upper surface, of width 2A, has infinite emissivity while 
the lower face and the ledges of widths A do not emit (see 
Fig.l). 2) A very strong strong magnetic field perpendic- 
ular to the electrodes inhibits the transversal components 
of electron velocities [6,8], but almost does not affect the 
total current [6,8,9]. 3) The emitted electrons leave the 
cathode with zero velocity [1,2,6]. 

If the cathode is in the {X, Z) plane and the mag- 
netic field in the F-direction the velocities v of electrons 
are parallel to the F-axis with mv'^{X, Y) ~ 2eU{X, Y), 
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where U{X,Y) represents the potential field while m, e 
are the electron mass and charge. The current density 
J{X), which clearly is K-independent, determines to- 
gether with v{X, Y) the density of electrons. Using the 
dimensionless variables 



X Y A ^ A 

X = y=— , a=— , A 



(1) 



2e y3/2 



J(A), 



the nonlinear Poisson equation for the potential then 
takes the form 



(2) 



The electron density p(x, y) and current j{x) are different 
from zero only in the shaded rectangle Q of Fig. 1 which 
shows a two dimensional cross section of our system. 
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FIG. 1. Geometry of the system 
Eq.(2), subject to the boundary conditions (BC), 
(p{x, 0) = for < a -|- A, 0(a;, 1) = 1 for \x\ < oo. 



(3) 



dy 



{x, +0) = for \x\ < a, 



is to be solved in the half-plane y < 1 to produce both 
functions 4'{x,y) and j{x) which are non- negative and 
symmetric about the y-axis. To do this we first solve 
eq.(2) approximately in the current region Q on a rather 
intuitive level. The problem is nonlinear here and it is 
not well posed if one disregards the field at > a. Con- 
sequently our solution will have a set of free parameters 
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which specify j{x) and (j}{x,y): in particular (j}{a,y) and 
|^(a~,j/). In the second step the potential (p{±a,y) is 
used as the BC and we obtain a Dirichlet problem for 
the Laplace eq.(2) in the outer region of the half-plane 
where j{x) = 0. We solve this problem using conformal 
mapping techniques and evaluate |^(a+,y). If one ex- 
cludes the points x = ±a, y = the electron density 
p{a^,y) is finite and p{a^,y) = 0, therefore the second 
derivative of (j){x, y) has a finite jump at x = a, while the 
first derivative must be continuous, i.e. 

^(«-,.) = |^«.^ri.o<v<l. (4) 

In the last step we satisfy approximately the matching 
condition (4) by adjusting the free parameters mentioned 
above using the least squares technique. This will give 
an approximate explicit form for j{x). 

The space charge region Q. We want to solve approx- 
imately eq.(2) where the function j{x) is not known 
nor are the BC for ^ at a; = ±a. When a = oo 
we have no x dependence and (2) becomes an ordi- 
nary equation which was solved in [1,2] yielding 4>i{y) = 
2/'*/^, ji{x) = 1. This gives the Child-Langmuir formula 
[1], Ji = (2e/m)i/2^3/2/g^£,2. For a » 1 it is reason- 
able to assume that j{x) ^ ji = 1 when a — |x| » 1 
and use also a stronger assumption that the difference 
(t)(x,y)/(j)\{y) — 1 is small almost everywhere (i.e. it does 
not exceed ~ 1 — 1.5 even near the edges of region). This 
difference however is not small at the cathode edges, 
x = ±a, where the electric field must match the field 
outside. The large gradients in the field lead to the ac- 
celeration of electrons and thus to a strong rise of the 
current density j{x) near the bomidary of the SCL flow. 

We represent (j){x,y) in the form y^^^[l + p{x,y)] and 
linearize the square root as [1 + /x(x, y)]~^/^ « 1 — 
fip{x, y), where the number /3 is chosen to minimize the 
integral of [l-/3,u-(l+^)^-^/^]^ on the interval < < 1. 
This yields /? ~ 0.328 with relative average error of ap- 
proximation around 2.2%. For fi = 0.2, 1, 1.5 the error 
is 2.36%.. 4.96%, 6.25%, 19.6% respectively We shall 
see later that for all A > 0.1 /x < 1.5. Substituting in (2) 
we obtain a linear equation 



y 



d'^H d^ii \ 8 ^1_ 
dx^ dy'^ J dy 
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(5) 



where we dropped a nonlinear term in the right side. The 
error due to this and to the linearization of the square 
root is negligible for small /z and decreases the right side 
by at most a factor ^ 0.7, in all the cases considered (see 
Table 1) including even /x ~ 2. 

Using the method of separation of variables we write 

^J'{x,y) = ^qifiix)ui{y), j{x) = l + j'^qifi{x), (6) 



with 



-ki{a—x) 



+ e 



-ki(a+x) 



\x\ < a. 



(7) 



Substituting (6) and (7) into (5) and assuming that 
(5), (3) are satisfied separately for each I = 1,2, ... gives a 
set of inhomogeneous equations 



, cPui 
dy'^ ' 3" dy 



y ^ + ^y— + {ky + )ui = i, (8) 



with the common BC ui{l) = 0. The parameters ki and 
qi will be determined later. The potential can be written 

in the form 



cl>{x,y)=y^/'+y^/'^qiMx)ui{y), 



(9) 



where the first term is the Child-Langmuir potential 
and the ui(y) are assumed finite. The relevant particu- 
lar solutions of (8), which can be expressed in terms of 
Lommel's functions S-i/6,i/(fciy), = ~ 16/3/6 [10], 
is given by the power series expansion 

2n 



«.(v) = j^E(-ir«.(^ 



(10) 



ao = 1, dn = 



0"n-l 



n2 + 5n/6+ (H-/3)/9 



As all ui{l) = the parameters ki are the increasing 
roots of (14): 3.881, 6.675, 10.065, 13.003, 16.316, 19.306, 
22.582, 25.600, 28.855, 31.891 for 1 < / < 10. They can 
be easily evaluated due to the rapid convergence of (10), 
asymptoticaly ki Iw. The free parameters g; will be 
used to satisfy (4). 

The 2-d mean current density over the whole cathode, 
which in terms of our scheme is given by 



i2=- rj(^)d^«^i+f El (1 

a Jo ki 



^ — 2kia 



), (11) 



is usually presented [8] as the 1-d current density ji = 1 
plus a correction: j2 = 1 + a /2a. Thus in the original 
units the mean current has the form 

J2 = Jl ( 1 + 1 , (12) 



W 



where W = 2A is the width of the cathode. Using (11) 
the parameter a is defined here by 



"=9E|;(i-^-^n 



(13) 



Electrostatic region. It seems clear that for a 1 the 
electric field in the vicinity of the boundary a; w a, < 

y<i 

is not affected much by the region x < —a, < y < 1, 
see Fig.l. This 

allows us to study a simpler electrostatic problem for 
a plane which is split according to Fig. 2(a). We modified 
a conformal transform in [11] to the form 

Z = 2'K~^[hl{y/w + Vw - I) - c^/w"^ - w], (14) 
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which maps the shaded half-plane z = x + iy on Fig. 2(a) 
onto the upper half-plane w = u + iv in Fig. 2(b). 
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FIG. 2(a). Plane z = x + iy FIG. 2(b). Plane w = u + iv 

Our Dirichlet problem with the BC on the real axis 
Qw = (which come from the previous section), 



{1, if -00 < u < 

^{a,y{u)), ifO<M<l 
0, ifu>l, 



(15), 



has the solution 



in the upper half-plane 
27r~^[arccos ^/u 



$(5,0) 



{u — s)^ + v"^ 



ds 



(16) 



Here by (14) y{u) 



c\Ju{l — u)]. The potential ^{u,v) ex- 
pressed in variables x,y represents (j){x,y) outside the 
space charge zone. Our task now is to match the inside 
electric field ff (a,y(w)) in the interval < u < 1 with 
the field outside 



dx 



[u, v = Q) 



TT ^u(i - u) a$ 

1 + c(l - 2u) 'dv 



(u,0). (17) 



Continuity of the electric field. The matching condi- 
tion (4) guarantees continuity of the electric field at the 
boundary between the space charge region Q with the 
vacuum. Using (7) and (9) we have at a; = a inside the 
space charge region the field intensity, 



dx 



^ = y4/3 



(18), 



which should be equal to the vacuum field (17). The 
exponentially small terms e"^*^'" can be dropped as a is 
assumed large. Both terms and ^ depend on all 
parameters qi , but in a different way. One cannot expect 
an exact equality because of the approximations made. 
We rewrite the matching condition (4) as 

G[y^^^] + E qi{G[y^^\i{y)] - hy^'My)} « 0, (19) 

where the functionals G can be written explicitly as in- 
tegrals with a logarithmic singularity. 

We minimize the least square divergence from zero of 
the expression (19) on the interval 0.15 < y < 0.85 with- 
out approaching the endpoints where our treatment is 
not entirely adequate. A standard procedure yields a set 
of linear algebraic equations for qi. We did not go further 
than Imax = 10. After the qi are computed one can find 
the current density (6) and the parameter a (13). 



The accuracy of this method can be evaluated to some 
degree by determining the relative average discrepancy 
A of electric fields at the boundary of the space charge 
region Q on the chosen interval of y. The results of our 
computations are shown in the Table 1, where for differ- 
ent values of the ledge A one can see also a(A) , parameters 
qi, fJ-max (near y = Q and x = ±a), and the relative height 
(see Fig. 3) of the current wings jmax at x = ±a. When 
we extend the interval of matching the electric fields up 
to (0.01,0.99) the quantities in the table stay approxi- 
mately the same, only A increases. This confirms the 
general validity of our method and simultaneously 

shows that the computation of electric fields near the 
corners of the rectangle Q is not very good. In particular, 
in the worst case (the most severe cathode regime, see 
also [6,7,12]) when A = 0, the electric field is singular at 
the cathode edges. The computation becomes unstable, 
we cannot therefore the data of Table 1 to be accurate 
there when the linearization fails too. 
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.0037 


.0032 
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.0142 
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TABLE 1 

When I runs from 1 to 10 the values of qi decrease 
approximately by a factor 20-50. Therefore a is evaluated 
very well by (13) where the kj increase from ~ 4 to 32. 
The accuracy of j{x) and fi{x) might be improved if one 
truncates (9) at a larger Imax though we would not expect 
dramatic changes. The cases A = oo (when the mapping 
is exact) and A = 1 are very close which means that 
the parameter a (a, A) as well as qi are approximately 
independent of a when a > I. Keeping the exponential 
terms in fj,{x, y) in the matching conditions (20) does not 
complicate the calculation and it will give 

only insignificant corrections. The current density j{x) 
at a; = increases in this case by about 2qle~^^^°^. If a 
is smaller, but 2A + 2a > 1, the scheme of computation 
is the same though the qi become functions of a and one 
cannot decrease a too much because the first term in (9) 
needs corrections. 

An important part of our analysis is the form (11) of 
fi{x) which implies that the current density (10) in a 
narrow region of width ~ 1 (£) in the original units) at 
the cathode edges has a sharp peak which decays faster 
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than cxp[— 3.88(a — Everywhere else j{x) is close 

to the 1-d current ji (x) = 1 with the exponentially small 
corrections. For illustration plots of the current density 
distribution (6) are shown in Fig. 3 for different widths of 
the cathode. 



FIG. 3. The current densities when A = 0.1, cl — cl, c2 — 
c2, c3 — c3 and c4 correspond to 2a = 8, 4, 2, 1 respectively. 

We can compare our curve c3 for a = 1 with PIC sim- 
ulations presented in [7]. There for zero cathode recess 
(dx = 0), ^ = 50 mm, and unfortunately unspecified 
width of the shroud a reasonable fit would be jmax = 3.9 
in our Table for A = 0.1 versus 3.2 in [7]. We get the 
half- width of the current density peak ~ 1.2 mm while 
in [7] it was 1 mm. Our magnetic field is stronger and 
we think also that PIC simulations with finer grids are 
closer to our computation, but diverge from experimental 
results because the real cathodes with their finite thick- 
ness and roundness do not have the very strong electric 
field intensities present however in the models. 

Generalization. We expect that this pattern of narrow 
wings of the current density holds also for finite flat cath- 
odes with perimeter P where the boundary region will 
have an area ~ PD if we assume reasonable restrictions 
on the curvature and return back to the original units. 
In the general case of a cathode with area S » PD the 
total current can be written as the sum I = SJi + PDJ. 
The " edge" current J, which is assumed here to be inde- 
pendent of geometry, can be viewed as distributed over 
the edges of the cathode of width equal to the distance 
between the electrodes. The ratio J/Ji can be evaluated 
in terms of the parameter a defined in (13). Comparing 
J = I/S= Ji(l + JPD/JiS) with eq.(12) for our geom- 
etry, where J = J2, S = 2AL and P = 2L (the length L 
of the cathode segment is arbitrary), we have J = J\a/2 
and finally 

^=^1(1 + "^), (20) 

which should be applicable in general situations. 

In particular the factor PD/2S in (20) becomes 
D/R for a circular cathode of the radius R and 
2DE{^Jl - C^/B"^) /ttC for an elliptical cathode with the 
half-axes B > C, where E{k) is the complete elliptical 
integral. For a rectangular cathode with the sides L and 
H it is equal to D{L-^ + H''^). 



Conclusions. 1) The current wings, Fig. 3, resemble 
simulated ones [6,7,12]. They arc high when the width of 
ledges A is small and the vacuum electric field near the 
cathode edges is strong. Their form becomes practicaly 
constant when the ledges are wider than the distance D 
between electrodes. 2) The shape of the current wings, 
which is determined by eigenvalues ki of (8), is roughly 
exponential and the 1-d current is restored up to a few 
percents at the distance D from the edges. 3) The param- 
eter a, which defines the net current density, depends on 
the width of ledges. An approximate empirical formula 



a(A) « 0.19 + 0.48e 



-3.7A 



(21) 
(For a 



agrees with the data in Table 1 within ~ 3.3%. 
different model with the constant current density a was 
estimated in [8] as close to 0.31.) 4) Our techniques of 
matching the electric fields at the boundary of the space 
charge region and using rather modest variations of the 
potential in the x direction is effective for approximate 
modelling the 2-d and 3-d flows of charged particles. 
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